NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

DEVELOPMENT OF GRAPHICAL POLE- ZERO, 

ROOT-LOCUS, BODE, NYQUIST , AND NICHOLS 
RESPONSES USING THE OPTSYSX PROGRAM 

’ay 

Michael Henry Laptas 
September 1 984 

Thesis Advisor: Daniel J. Collins 

Approved for public release; distribution unlimited 

T222961 




SECURITY CLASSIFICATION of THIS ^AGE Data Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subtitle) 

Development of Graphical Pole-Zero, 

Root-Locus, Bode, Nyquist, and Nichols 
Responses Using the OPTSYSX Program 


S. TYPE OF REPORT 6 PERIOD COVERED 

Master 1 s Thesis 
September 1 984 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR)'*.) 

Michael Henry Laptas 


8, CONTRACT OR GRANT NUMBERfaJ 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93943 


10. PROGRAM ELEMENT. PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey , California 93943 


12. REPORT DATE 

September 1 984 


13. NUMBER OF PAGES 

21 9 


14. MONITORING AGENCY NAME A ADDRESS^// different from Controiiing Office) 


IS. SECURITY CLASS, (ot this report) 


IS*. DECL ASSI F| CATION/ DOWNGRADING 

schedule { 


16. DISTRIBUTION STATEMENT (ot this Report) 

Approved for public release; distribution unlimited 


17. DISTRIBUTION STATEMENT (of the abstract entered in Block 20, if different from Report) 


18. supplementary NOTES 


19. KEY WORDS (Continue on reverse aide if necessary and identity by block number) 

OPTSYSX Nyquist Noise 

Pole-Zero Nichols Compensator 

Root- Locus Transfer function 

Bode Open Loop 


20. ABSTRACT (Continue on reverse aide It necessary and identify by biock number) 

This thesis discusses the modification of and additions 
to an existing Optimal Systems Control FORTRAN Program (0PTSYS) 
originally obtained from Professor Arthur E. Bryson of Stanford 
University. This program has been subsequently redesigned to 
run interactively on the IBM 3033 VM/CMS by Lieutenant Commander 
John G. Hoden, and additions by Commander Harry A. Diel provide 
the user with a highly accurate graphic time response to a 
system designed using the OPTSYSX program. 



DD , jan^73 1473 



EDITION OF 1 NOV 6S IS OBSOLETE 

S N 0102- LF- 01 4- 6601 



SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 



SECURITY CLASSIFICATION OF THIS PAGE (1Th*n Enlatad) 



The addition of the FORTRAN program OPTGRAPH gives the 
user the capability to obtain classical analysis (Pole-Zero Map, 
Root Locus, Bode, Nyquist, and Nichols) of transfer functions, 
calculated for a system designed with the OPTSYSX program. The 
OPTGRAPH program uses high resolution precision plotting 
software to provide the user with highly accurate frequency 
response plots. 

This series of programs permits users to rapidly carry out 
simulation, analysis, and design of Optimal Systems Control 
problems in a totally interactive mode. 



S'N 0102- LF- 014- 6601 



2 

SECURITY CLASSIFICATION OF THIS PAGE(W*n Data Eniarad) 



Approved for public release; distribution unlimited 



Development of Graphical Pole-Zero, Root-Locus, Bode 
Nyquist, ana Nichols Responses using the 
OPTSYSX Program 



by 



Michael H. Laptas 
Lieutenant, United States Navy 
B. S. , Purdue University, Lafayette, 1977 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 



from the 



NAVAL POSTGRADUATE SCHOOL 
September 1984 



™ UBRARY 

5fe TG . ■ D!JAT E SCHOOL 
1 • ^UFORNIA 93943 



ABSTRACT 



This thesis discusses the modification of and additions 
to an existing Optimal Systems Control FORTRAN Program 
(OPTSYS) originally obtained from Professor Arthur E. Bryson 
of Stanford University. This program has been subsequently 
redesigned to run interactively on the IBM 3033 VM/CMS by 
lieutenant Commander John G. Hoden, and additions by 
Commander Harry A. Diel provide the user with a highly accu- 
rate graphic time response to a system designed using the 
OPT S Y SX program. 

The addition of the FORTRAN program OPTGRAPH gives the 
user the capability tc obtain classical analysis (Pole-Zero 
Map, Root Locus, Bode, Nyguist, and Nichols) of transfer 
functions, calculated for a system designed with the OPTSYSX 
program. The OPTGRAPH program uses high resolution precision 
plotting software to provide the user with highly accurate 
frequency response plots. 

This series of programs permits users to rapidly carry 
out simulation, analysis, and design of Optimal Systems 
Control problems in a totally interactive mode. 



4 



ouoley. knox ^ School 

* .•.’. , ' 9S343 

TABLE OF CONTENTS 



f/iTS t ^ 

C • / 



I. INTRODUCTION 9 

II. THE OPTSYSX COMPUTER PROGRAM 11 

A. BACKGROUND 11 

B. OVERVIEW 11 

C. OPTSYSX MODIFICATIONS 12 

D. OPTSYS EXEC MODIFICATIONS 13 

III. THE OPGRAPH PROGRAM 15 

A. PROGRAM OVERVIEW 15 

1. Program Language 15 

2. Graphics Package 15 

3. Library Subroutines 15 

4. Program Composition 16 

B. GENERAL PROGRAM OPERATION 16 

1. User Protection Features 16 

2. Large Order System 17 

C. SYSTEM/MODEL DESCRIPTION 18 

D. PROGRAM ORGANIZATION 19 

1. Main Program 19 

2. Program Sequencing 20 

3. Graphic and/or Tabular System Response . . 21 

E. INTERACTIVE EXAMPLES 23 

F. CHARACTERISTIC EQUATION ANALYSIS 24 

1. Pole-Zero Map 24 

2. Pole-Zero Map (Interactive Example) ... 25 

3. Root-Locus 33 

4. Root-Iccus (Interactive Example) 34 



5 



G. FREQUENCY RESPONSE ANALYSIS 43 

1. Bode Response 44 

2. Bode Response (Interactive Example) ... 45 

3. Nyquist Response 51 

4. Nyquist Response (Interactive Example) . . 52 

5. Nichols Response 52 

6. Nichols Response (Interactive Example) . . 52 

IV. CONCLUSIONS AND RECOMMENDATIONS 59 

A. CONCLUSIONS 59 

B. RECOMMENDATIONS 60 

1. Root Finding Routines 60 

2. Program Memory Requirements 61 

3. Further Modifications 61 

APPENDIX A: OPTSYS EXEC LISTING 62 

APPENDIX B: OPTSYSX PROGRAM LISTING 71 

APPENDIX C: OPTGRAPH PROGRAM LISTING 160 

APPENDIX D: OPGRAPH LISTING 215 

LIST CF REFERENCES 217 

EIBL IOGRAPHY 218 

INITIAL DISTRIBUTION LIST 219 



6 



LIST OF FIGURES 



3. 1 Pole-Zero Map Example 53 

3.2 Root Locus Response Example 54 

3.3 Bode Response Example (Magnitude) 55 

3.4 Bode Response Example (Phase) 56 

3.5 Nyguist Response Example 57 

3.6 Nichols Response Example 58 



7 



SYMBOLS 



A 

B 

C 

D 

F 

G 

GAM 

H 

K 

Nc 

Ng 

Ns 

No 

Q 

R 

c 

U 

V 

w 

wO 

x 

xdot 

xe 

xedot 

y 

2 



State (Ns, Ns) or Output (No, No) Weighting Matrix 
Control (NC/Nc) Weighting Matrix 
Control Gain Matrix (Nc,Ns) 

Control (No,Nc) or Noise (No,Ng) Feedforward 
Matrix 

Open-Loop Dynamics Matrix (Ns, Ns) 

Control Distribution Matrix (Ns,Nc) 

State Disturbance Distribution Matrix ( N s , N g ) 
Measurement Scaling Matrix (No, Ns) 

Estimator Gain Matrix (Ns, Nob) 

Number of Controls 

Number of Process Noise Sources 

Number of States 

Number of Observations or Measurements 
White Process Noise Covariance Matrix (Ng,Ng) 
White Meas. Noise Covariance Matrix (No, No) 
Steady-State Covariance Matrix of Control (Nc,Nc) 
Control Vector (Nc,1) 

White Measurement Noise Vector (No,1), with 
Zero Mean and Covariance Matrix R 
White Process Noise Vector (Ng,1), with 
Zero Mean and Covariance Matrix Q 
Constant Disturbance Vector (Ng, 1) 

State Vector (Ns,1) 

Derivative cf State Vector (Ns,1) 

Estimate of State Vector (Ns, 1) 

Derivative cf Estimate of State Vector (Ns,1) 
Output Vector (No, 1) 

Measurement Vector (No, 1) 
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I. INTRODUCTION 

The purpose of this thesis is to describe the modifica- 
tion and additions tc the existing FORTRAN program (OPTSYS) 
which is used in the study and application of Optimal 
Systems Control theory. 

The optimal SYstems control program was originally 
developed by Hall [Ref. 1 ] to support his research in 
rotary-wing aircraft control systems. Later program modifi- 
cations were made by Walker [Ref. 2] and Liu [Ref. 3] of 
Stanford University, and are designated OPTSYS 4 and OPTSYS 
5 respectively, OPTSYS modifications made by Hoden [Ref. 4] 
were primarily devoted to creating a user-friendly interac- 
tive version (OPTSYSX) of the OPTSYS 4 program. The latest 
modifications by Diel [Ref. 5] allowed the user to save 
matrices for subsequent runs of the OPTSYSX program, and 
formed data file sets for the time response program OPTCALC. 

The intent of this thesis work was to develop an inter- 
active program to plot the Pole-Zero map, Root-Locus, Bode, 
Nyquist, and Nichols responses to the open loop, closed loop 
noise, and compensator transfer functions of a State 
Variable Control System which has been developed using 
OPTSYS Program. Miner modifications to the OPTSYSX Program 
were necessary to construct a data file sets for plotting 
the Pole-Zero Map, and in calculations for the Root-Locus, 
Bode, Nyguist, and Nichols plots. 

It is assumed that the user is familiar with the basic 
concepts cf Control Theory and Optimal Systems Design. The 
symbol/naming conventions of Bryson [Ref. 6]. are used in 
the discussion of program operations and descriptions of 
problems using the OPTSYS System. 
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An overview of the OPTSYSX Program capabilities and of 
modifications to the existing programs is presented first. 
This is followed by a description of the program (OPTGRAPH) , 
which was developed tc plot the Pole-Zero Map and to perform 
calculations for the Root-Locus, Bode, Nyguist, and Nichols 
plot s . 

This work concludes with examples of various types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final results. 
Complete program listings for the OPTSYS EXEC program, 
CPTSYSX program, and CPTGEAPH program are included in 
appendices A,B, and C respectively. 
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II. THE OP TSYS I COM PD TER PROGRAM 



A. BACKGROUND 

OPTSYSX is a double-precision/ interactive FORTRAN 
program employing modern control theory analysis techniques. 
Its capabilities include the calculation of the cjen-loop 
eigensystem, and the stationary closed loop system; the 
synthesis of regulators and filters; along with power spec- 
tral density, and modal distribution computations. The 
modifications introduced to the OPTSYSX program by this 
thesis work do net affect the program's original 

capabilities. 



B. OVERVIEW 

OPTSYSX is an extremely large and complex program with 
over 3000 lines of code. To use the program in its small 
version (dimensioned for a 32 X 32 "F" , "G”, and "H" 

matrices), the user must extend his Virtual Machine (VM) 
memory capacity beyond the default VM memory capacity of 720 
kilobytes to 1024 kilobytes. 

Any significant increase in the OPTSYSX program size, 
with the the resulting requirement for additional memory 
capacity, would cause user difficulties. For this reason 
the task of obtaining the Pole-Zero maps, Root-Locus, and 
Bode, Nyquist, Nichols plots for the open loop, noise and 
compensator transfer functions was relegated to a separate 
program (OPTGRAPH) . Three data files (OPTGROL DATA, OPTGRNO 
DATA, and OPTGRCM DATA) , containing the open loop, noise, 
and compensator transfer functions, are formed by the OPTSYS 
program to be transfered to the OPTGRAPH program. The 
OPTSYS EXEC program provides an interface between the 



OPTS Y SX program, the time response programs (OPTCALC and 
OPTPICT) , and the OPTGRAPH program. 

C. OPTSISX MODIFICATIONS 

CPTSYSX program modifications consist of the addition of 
write statements, three flags, and a short routine to inter- 
pret user inputs. 

Write statements, to input data to the OPTGROL, CPTGRNO, 
and CPTGRCM DATA files were added to the main program, 
subroutine INNER, subroutine ZEROS, and subroutine POLES. 
System information consisting of number of states (Ns) , 
number of controls (Nc) , number of measurements (No) , number 
of process noise sources (Ng) , type transfer function 
(ITFX) , Markov Parameter (IE) , and two flags to be used by 
the CETSYS EXEC and the OPTGRAPH programs is inputed from 
the main program. The write statements added to the subrou- 
tine INNER, and subroutine POLES input the poles for the 
open loop and noise transfer functions, and compensator 
transfer function respectively. The zeros, the numerator 
order, the gain, and the input and output numbers for all 
three transfer functions are obtained from the subroutine 
ZEROS. 

The first of the two flags, input with the system infor- 
mation from the main program to the three data files, serves 
as a marker for the CPTSYS EXEC to locate the correct data 
line for reading system inputs. The second flag, which is 
set by the routine tc interpret user inputs, serves as a 
signal for the OPTSYS EXEC to either load the OPTGRAPH 
program if any one of the three transfer functions was 
calculated by the OETSYSX program, and abort loading the 
OPTGRAPH program if none of the three transfer functions 
were calculated. The third flag (ITFX) is an existing 
OPTSYSX program flag. The flag (ITFX) is passed from the 
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subroutine TF to the subroutine ZEROS to identify the calcu- 
lation cf either an open loop, noise, or compensator 
transfer function zeros. 

The short routine added to the main program of OPTSYSX 
reads the flags ITF1, ITF2, ITF3, and IRET. The flag sent 
to the CPTGROL , OPTGRNO, or OPTGRCM DATA files is set to 
allow loading of transfer function data file to the CPTGRAPH 
program if the user selects either of the transfer function 
options number two (calculate poles, residues, and zeros - 
ITF 1 = 1, ITF2 = 1, or ITF3 = 1) number three (calculate 
poles, and zeros - ITF1 = 2, ITF2 =2, or, ITF3 = 2) when 
calculating the open loop, noise, or compensator transfer 
functions. The selection of any other transfer function 
option will not calculate the transfer function zeros. The 
flag IRET is read to determine if the user desires to rerun 
a problem without exiting the OPTSYSX program. Should the 
user elect to rerun a problem without exiting OPTSYSX (IRET 
= 1) , the three data files are set back to the beginning 
(REWIND) to accept the new problem system data. 

D. OPTSYS EXEC MODIFICATIONS 

The OPTSYS EXEC was written by Diel [Ref. 5] to inter- 
face between the OPTSYSX program and the time response 
programs 0PTCA1C and CPTP10T. The modifications to the exec 
program were to add the option to select running the 
OPTGRAPH program to the exec menues and a short section to 
interface the OPTGRAPH program with the OPTSYSX program. 

The added section to interface between the OPTSYSX and 
OPTGRAPH programs defined loader size (LDRTBLE) , file defi- 
nitions (FILEDEF) , and text libraries (TXTLIB) for the 
CPTGRAPH program. The first line of the three data files is 
read to ensure that there has been at least one of the 
three transfer functions calculated by the OPTSYSX program 
prior to loading OPTGRAPH program. 
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The OPTSYS EXEC routine is also used to send plots to 
the VEESATEC plotter from the OPTPLOT program and the 
OPTGEAPH program. 
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III. THE OPGRAPH PROGRAM 
A. PROGRAM OVERVIEW 

CPTGRAPH is an interactive FORTRAN program which uses 
the transfer function poles and zeros calculated by OPTSYSX 
to calculate, tabulate, and plot the Pole-Zero Map, Pode, 
Nyguist, Nichols, and Root-Locus responses for either the 
open loop, closed loop noise, or compensator transfer 
functions. 

1 • Pro gra m Lang uage 

CPTGRAPH is programmed in FORTRAN following the 
conventions of the IBM System /360/370/ FORTRAN IV language. 
CPTGRAPH has been compiled and run under both FORTRAN IV 
(G 1 ) and FORTRAN H (extended) compilers on the IBM 3033. 

2 • G raphics Pack age 

This program uses the Display Integrated Software 
System and Plotting System Software (DISSPLA) developed and 
distributed by the Integrated Software Systems Corporation 
(ISSCC) of San Diego, California. The DISSPLA package is a 
library cf FORTRAN subroutines using FORTRAN IV conventions. 

3 • lib rary Sub r ou tines 

CPTGRAPH uses the International Mathematical and 
Statistical Library (IMSL) subroutine ZRPOLY to calculate 
the roots of a polynomial equation. ZRPOLY is a double 
precision FORTRAN subroutine to which is input a polynomial 
with real coefficients (double precision) in terms of 
decreasing powers of the polynomial. The subroutine finds 
the rccts of the polynomial and returns a double precision 
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complex array containing the roots. ZRPOLY is capable of 
accepting polynomials greater than 0 order and less than 101 
order . 

4 . Pro gra m Co m position 

CPTGRAPH has one main program and 24 subroutines. 
The main program and its subroutines may be divided into 
four basic categories: 

1) File Data Input 

2) Interactive Data Input 

3) Plot Setup and Sequencing 

4) Calculation 

A brief and general description of the program and 
its subroutines will follow in subsequent sections. 

B. GENERAL PROGRAM OPERATION 

OPTGEAPH was written to satisfy two specific objectives. 
The first is to create a program to be used as a instruc- 
tional tool for students taking controls related courses, 
and the second is as method to assist in the analysis of 
actual large order control system problems. 

To accomplish these objectives the emphasis while 
writing this program was to make the program as user 
oriented as possible by attempting to eliminate ambiguities 
and providing features which would protect the user from 
inadvertent wrong entries. An effort was also made to mini- 
mize the Virtual Machine (VM) memory required to execute the 
program to facilitate the analysis of a large order system. 

1 * Oser Protection Features 

Three basic methods were used to protect the user 
from an abnormal program termination in the case of an inad- 
vertent wrong entry. The first was to construct the 
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subroutines to display a summery of the user inputs and 
allowing changes prior to leaving the subroutine. For the 
case where the user input is expected to fall within a given 
range, the input is examined and the user is either allowed 
to continue for an input within the range or issued an error 
or warning and reasked the question for an input outside the 
range. The final method is used by the three subroutines 
which read the user input from the screen. The subroutines 
RDINT, RDREAL,and RECHAR expect an integer input, real 
number input, and a logical "yes" or "no" respectively. In 
the event that a "null" line is inadvertently entered, the 
user is issued a warning and allowed another opportunity to 
enter the the correct input before abnormally terminating 
the program. The entry of two "null" lines also offers the 
user the option to to exit the OPTGRAPH program at ether 
than normal program exit points. The subroutine RDCHAR also 
examines the input and issues a warning if the input is not 
either "yes" or "no". 

2 • large Order S ystem 

The large order control problem being considered for 
analysis by the OPTSXSX and the OPTGRAPH programs is the 
X-29A aircraft longitudinal axis backup mode system. This 
system has a {98 X 98) »F" matrix , a (2 X 98) "H" matrix, 
and a {98 X 1) "G" matrix. 

To accommodate this large order system the OPTGRAPH 
program was dimensioned to accept a system with a maximum of 
99 states (Ns) , 12 controls (Nc) , 12 process noise sources 
(Ng) , and 12 measurements (No). The Virtual Machine (VM) 
memory requirements for the program were minimized by the 
reuse of memory locations allocated for storage of calcu- 
lated data after that data has been either tabulated or 
plotted. This method of memory reuse has allowed the 
OPTGRAPH program to be dimensioned to accommodate a large 
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order system and still operate with less than the 1024 
KILOBYTES VM memory required to run OPTSYSX in its small 
vers ion . 

C. SYSTEM/MODEL DESCBIPTION 

The system equations used by the OPTSYSX program for the 
calculation of the transfer functions are of the state vari- 
able form. The system equations are: 
system model 

xdot = [F]*x + [G]*u + [GAM ]*w (3.1) 

measurement equation 

z = [ H ]*x + [ D ]*u + v (3.2) 

estimator equation 

xedot = [ F ]*xe + [G]*u + [K]*(z - [H]*xe) (3.3) 

open loop transfer function 

[ H ]* | s[ I ] - [F]|-**[G] (3.4) 

closed loop noise transfer function 

[ H ]* | s[ I ] - [ F ] 1 - 1 *[ Gam ] (3.5) 

compensator transfer function from measurement to input 

[C]*(]s[I]-[F] + [G]*[C]+[K ]*[H ]|-i) *[K] (3.6) 



where 

u = Control vector (Ns X 1) 
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w = White process noise vector (Ng X 1) 

x = State vector (Ns X 1) 

xdot = Derivative of the state vector (NS X 1) 
y = Output vector (No X 1) 

z = measurement vector (No X 1) 

[C] = control gain matrix (Nc,Ns) 

[D] = control feed-forward distribution matrix (Nc, Nc) 
of the control vector (u) 

[F] = open-loop dynamics matrix (Ns, Ns) 

(System matrix or Plant matrix) 

[G] = control distribution matrix (Ns,Nc) 

[GAM] = state disturbance distribution matrix (Ns,Ng) 

[H] = measurement distribution matrix (No, Ns) 

[I] = identity matrix (Ns, Ns) 



D- PBOGBAM ORGANIZATION 
1 • Mai n Pr ogra m 

The main program presents a menu allowing the user 
to select for analysis either the open loop transfer func- 
tion, the closed loop noise transfer function, the compen- 
sator transfer function, or exit the OPTGRAPH program. The 
data file for the selected transfer function is examined by 
the main program to ensure that the file contains the 
required transfer function data. Should the data file be 
incomplete the program gives the user a warning message 
indicating the nature of the missing data and allows the 
user tc select another option. 

The main program menu: 

OPTGRAPH 

DO 100 DESIRE TO AHALYZE: 
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1. OPEN LOOE TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER OPTION NUMBER. 

Data is read from the selected file by the main 
program and formed into three arrays and seven integer vari- 
ables. The exchange of data between the main program and 
the subroutines, and between the subroutines is by the 
extensive use of labeled common statements. The data trans- 
fered through the subroutine call statements is limited 
primarily to flags, and constants. The program moves to its 
primary sequencing subroutine (GRAPH) upon the completion of 
file data entry. 

2 . P rogra m Seq uencing 

The subroutine GRAPH first presents the user with a 
summary cf the the transfer system received from the data 
file. The user is then asked which transfer function he 
desires to to analyze using the conventions of the OPTSYSX 
program describing the transfer functions by input number 
and output number. A menu is presented which allows the 
user to select graphical and/or tabular system response 
(Pole-Zero map. Root-locus, Bode, Nyguist, or Nichols) , 

select another transfer function, or exit the OPTGRAPH 

program. Upon completion of a graphical and/or tabular 
response option, or change of transfer function option, the 
program returns to the GRAPH menu to allow the user to 
choose another option or exit to the main program. 

Subroutine GRAPH data summary and menu: 

OPTGRAPH 

THE OPEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS CONTAINS: 
1. HUMBER OF STATES = 4 
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2. NUMBER OF CONTROLS (INPUTS) = 1 

3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 
H. MARKOV PARAMETER = 10.0**- 6 

CLEAR SCREEN TO CONTINUE 



OPTGRAEH 

AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, AND/ 
OR NICHOIS PLOT IS DESIRED FOR: 

INPUT # ? 

OUTPUT # ? 



OPTGRAEH 

OPEN LOOP TRANSFER FUNCTION 
INPUT # = 1 

OUTPUT f = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR DATA FOR: 

1. POLE-ZERO MAP 

2. ROOT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPU1/OUTPDT COMBINATION (SAME TRANSFER 
FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT OPTGRAPH 
ENTER OPTION NUMBER. 

3. Gra p hi c and/or Tabular System Response 

The selection of a graphical and/or tabular system 
response option for an Pole-Zero map, Root-Locus, Bode, open 
loop Nyguist, or open loop Nichols calls the subroutines 
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PZERO, RT10, BODE, NYQST, or NICHOL respectively. While the 
response of these subroutines is different, their organiza- 
tion is similar. 

The subroutines PZERO, RT10, BODE, NYQST , and NICHOL 
primarily act as sequencing and plot setup subroutines. The 
interactive questioning of the user for graphic/tabular 
response information and the system response calculation has 
been relegated to other subroutines. The subroutines PZERO, 

. RTLO, ECDE , NYQST, and NICHOL are divided into two major 
sections with the first being devoted to producing a graph- 
ical output and the second to producing a tabular output. 

The user is presented a subroutine menu offering the 
option for graphic response, tabular response, or exiting 
the subroutine. The program moves to either the first 
section (graphic response) , second section (tabular 

response) or back to the subroutine GRAPH menu depending on 
the option selection. The program returns to the subroutine 
menu upon completion cf the of either the graphic response, 
or tabular response section, allowing the user to select 
another response option, or exit the subroutine. 

Features that are common to the five subroutines 
will be described here and features that are unique to one 
of the subroutines will be described in the following 
section under the appropriate response heading. 

For a graphic response, the user is given the choice 
of two printer options (TEK 618, or VERSATEC) . With the 
selection of the TEK 618, the plot page size is defaulted to 

11 inches by 8.5 inches, and the plot will be presented on 

the TEK 618 screen. The VERSATEC option gives the user the 

option to create up to a 21 inch by 21 inch plot (maximum 

for the VERSATEC printer) and causes a DISSPLA METAFILE to 
be created. A scaling subroutine (PSCALE) scales the plot 
heading, legend, etc. in proportion to the page width 
selected. The VERSATEC option offers the the advantage of a 
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high quality print from the VERSATEC printer 
VERSATEC printer option in the OPTSYS EXEC. The 
a DISSP1A METAFILE also gives the user (aft 
optsys) the opportunity to use printers other 
618 or VERSATEC. A major disadvantage with 
option is that the user must exit the OPTGR 
fcefore he can print a graphic system response, 
reason, it is recommended that the system respo 
plotted on the TEK 618 prior to selecting t 
option. 



using the 
creation of 
er he exits 
than the TEK 
the VERSATEC 
APH program 
For this 
nse first be 
he VERSATEC 



For a tabula 
which gives him t 
printer, or disk) t 
asked if he desires 
data has been calcul 
is presented a menu 
the tabular data rou 
Examples of 
following sections 
system responses. 



r response, the user is presented a menu 
he choice of three devices (screen, 
o send the tabular output. The user is 
to make any changes, after the tabular 
ated. If the answer is affirmative, he 
and allowed to make changes and rerun 
tine. 

these features are included in the 
with the interactive examples for the 



E. INTERACTIVE EXABPIES 

The open-loop, closed loop noise, and compensator 
transfer functions are are handled identically by the 
OPTGR AP H program, with the only difference being the program 
and graph headings identifying the transfer function type. 
The frequency response programs for the Bode, Nyquist, and 
Nichols responses use identical setup and tabular data 
routines. 

A good overall representation of the OPTGRAPH program 
capabilities will be demonstrated by single interactive 
terminal session examples for the pole-zero, root locus, and 
bode responses using selected transfer functions. Tabular 
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data examples will te provided for the Pole- Zero, Root 
Locus, and Bode responses, and graphic response examples 
will te provided the Pole-Zero, Root Locus, Bode, Nyguist, 
and Nichols responses. 

F. CHARACTERISTIC EQUATION ANALYSIS 

The subroutines PZERO and RTLO are similar in their 
basic formats. The pole and zero locations for the selected 
transfer function input and output numbers are read into two 
single dimension arrays. The extraneous zeros, calculated 
by the OPTSYSX program, are eliminated by comparing the 
zeros to the Markov parameter sent with the system informa- 
tion from the OPTSYSX program. For a graphical response, 
the user is asked to define the plot limits in terms of 
X-coordinates (real axis) and Y-coordinates (imaginary 
axis) . Points for the plot which fall outside these limits 
will be ignored. System information (tranfer function 
input number, output number, and (DC) gain) is listed in the 
graph legend. 

1 . Po l e-Z ero Ma p 

The subroutine PZERO converts the double precision 
numbers for pole and zero locations to single precision for 
compatibility with the with the Graphics Package and plots 
the poles and zeros within the plot limits defined by the 
user for the graphical response routine. For tabular data 
response the user is given system information (numerator 
order, denominator order, and transfer function (DC) gain) 
for the selected transfer function input number and output 
number, and the pole and zero locations. 
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2 • Fo i e-Z e ro Map ( Inte r activ e Examp le ) 

The following example of a simplified open loop transfer 
function for an airplane with an autopilot in the longitu- 
dinal mode was obtained from [Ref. 7,p.315]. 

The aircraft system is represented by: 

[D] = [0] 



0. 


i. 


• 

o 


0. 1 


• 

o 


0. 


1. 


0. 1 


V 

o 


0. 


0. 


1. 1 


0. 


16. 


-12. 


-3. | 



I 0. | 

[G] = | 0. | 

I 1. I 

1 - 2 . | 



[GAM ] = [0 ] 



[ H ] = [ 1. 0. 0. 0. ] 



a. Example cf Graphic Response (Terminal Session) 

BEGIN RECORDING OF TERMINAL SESSION 
R; T=0.0 1/0.02 19: 40: 17 
EXECUTION BEGINS 



OPT GRAPH 

DO YOU DESIRE TO ANALYZE: 
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1. OPEN LOOP TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER OPTION NUMBER. 

1 



OPTGRAPH 

THE OPEN IOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS 
CONTAINS: 

1. NUMBER OF STATES = 4 

2. NUMBER OF CONTROLS (INPUTS) = 1 

3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 

4. MARKOV PARAMETER = 10.0**- 6 

CLEAR SCREEN TO CONTINUE 
OPTGRAPH 

AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST , 
AND? OR NICHOLS PLOT IS DESIRED FOR: 

INPCT # ? 

1 

OUTPUT # ? 

1 



OPTGRAPH 

OPEN LOOP TRANSFER FUNCTION 
INPCT # = 1 

OUTPUT # = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TAEULAR 
DATA FOR: 

1. POLE-ZERO MAP 
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2. ROOT-LOCU S 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. KICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 

ENTER OPTION NUMBER. 

1 



POLE-ZERO MAP 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER OPTION NUMBER. 

1 



POLE-ZERO MAP 

OPEN LOOP TRANSFER FUNCTION 
PLOTTER SELECTION AND PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SE LECTION OF TEK618) 

WHICH PLOTTER DO YOU DESIRE: 

1. TEK618 

2. VERSATEC 
CHOOSE OPTION 1 OR 2 

2 

PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 

HEIGHT = 

6 
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WIDTH 



8 

PLOTTING LIMITS FOE GRAPH 
X AXIS (REAL AXIS) 

X MINIMUM = ? 

-15 

X MAXIMUM = ? 

5 

Y AXIS (IMAGINARY AXIS) 

Y MINIMUM = ? 

-5 

Y MAXIMUM = ? 

5 

PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 

X MINIMUM = -15.00 

X MAXIMUM = 5.00 

Y AXIS (IMAGINARY AXIS) 

Y MINIMUM = -5.00 

Y MAXIMUM = 5.00 

DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES” OR "NO", 
n 

PLOT HEADING 

HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3 MAX ) 

ENTER NUMBER OF LINES. 

(ENTER "0" FOR NO HEADING) 

3 

A MAXIMUM OF 32 CHARACTERS PER LINE IS 

ALLOWED 

LINE 1 IS: 

aircraft with an autopilot in 
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LINE 2 IS 



the longitudinal mode 
LINE 3 IS: 
ref. 7; page 315 

PLOT HEALING 

LINE 1 : AIRCRAFT WITH AN AUTOPILOT IN 

LINE 2 : THE LONGITUDINAL MODE 

LINE 3 : REF. 7; PAGE 315 

DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES" OR "NO", 
n 

>USI NG A PRE- ALLOCATED DATASET FOR UNIT FT17F001 
>USI NG A PRE- ALLOCATED DATASET FOR UNIT FT18F001 



A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 



POLE-ZERO MAP 

OPEN LOOP TRANSFER FUNCTION 

DO YOU DESIRE TO MAKE ANY CHANGES TO: 

1. PLOTTER / PAGE SIZE 

2. GRAPH LIMITS 

3. HEADING 

4. NO ADDITIONAL CHANGES - PLOT POLE-ZERO 
MAP 

5. NO CHANGES - EXIT POLE-ZERO PLOTTING 
ROUTINE 

ENTER OPTION NUMBER. 

5 



POLE-ZERO MAP 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 
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1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER OPTION NUMBER. 

3 



OPTGRAPH 

OPEN LOOP TRANSFER FUNCTION 

INPUT # = 1 

OUTPUT # = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. POLE-ZERO MAP 

2. ROOT-LOCUS 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
OPTGRAPH 

ENTER OPTION NUMBER. 

7 



OPTGRAPH 

DO YOU DESIRE TO ANALYZE: 

1. OPEN LOOP TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER OPTION NUMBER. 

4 
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E; T=2 . 8 0/4. 67 1 9: 45 : 15 
record off 

END RECORDING OF TERMINAL SESSION 

The preceding example follows at the end of the 
chapter as figure 3. 1 . 

b. Example of Tabular Data (Terminal Session) 

The entering and exiting option menus for the 
tabular data terminal session are identical to the graphic 
response terminal session option menus, and have been elimi- 
nated from the following terminal session. 

BEGIN RECORDING OF TERMINAL SESSION 
R; 1=0.01/0.02 19:51:03 
EXECUTION BEGINS. . . 

POLE-ZERO MAP 

OPEN ICCP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER OPTION NUMBER. 

2 



POLE-ZERO MAP 

OPEN LOCP TRANSFER FUNCTION 

DO YOU DESIRE TABULAR OUTPUT TO GO TO: 

1. SCREEN 

2. PRINTER 

3. DISK (OPGRAPH LISTING) 

ENTER CPTICN NUMBER. 

1 
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POLE-ZERO HAP 

OPEN LOOP TRANSFER F DNCTION 
INPUT NUMBER = 1 

OUTPUT NUMBER = 1 



TRANSFER FUNCTION 

DENOMINATOR ORDER 
POLS LOCATIONS 

REAL PART 

0.0000D+00 
0.1000D+01 
-0.2000D+01 
-0.20 00D + 0 1 

NUMERATOR ORDER = 
ZERO LOCATIONS 

REAL PART 
-0.1000D+01 



(DC) GAIN = 0.1 000D+01 

= 4 

IMAGINARY PART 
0. 0000D+00 
0.0000D+00 
0.3464D+01 
-0.3464D+01 

1 

IMAGINARY PART 
0.0000D+00 



DO YOU DESIRE TO CHANGE OUTPUT DEVICE? 

TYPE "YES" OR "NO", 
n 



POLE-ZERO MAP 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 

ENTER OPTION NUMBER. 

3 

R; T=1. 60/3.23 19:53:00 
record off 



32 



END RECORDING OF TERMINAL SESSION 



3 . Root-Loc us 



The subroutine RTLO creates a closed loop system 
(with negative unity feed-back) from the selected transfer 
function input numbers and output numbers. Transfer func- 
tion poles and zeros are formed into two separate polyno- 
mials by the subroutine MAKPOL. The numerator polynomial 
(zeros) is multiplied by the gains, and the two polynomials 
are added to form a single polynomial. Complex roots of the 
of the polynomial are calculated by the IMSL library subrou- 
tine ZRPCLY. The complex roots are separated into real and 
imaginary components for either plotting or tabulating. 

The plotting routine first plots the transfer func- 
tion pole and zero locations then plots the root locations 
as they are received from the subroutine ZRPOLY. For the 
plotting routine the gain interval specified by the user is 
divided into 2000 evenly spaced points for the subroutine 
ZRPOIY to calculate roots. The poles, zeros, and roots are 
plotted within the plot limits specified by the user. 

For the tabular data routine the user is given the 
option tc choose between 1 and 500 points for the subroutine 
ZRPOLY to calculate roots. The tabular data output provided 
for the user is system information (numerator order, denomi- 
nator order, and transfer function (DC) gain) for the 
selected transfer function input number and output number, 
the pcle and zero locations, and the root locus gains and 
roots 



A positive unity feedback system may be 
specifying a negative gain range vice a positive 



analyzed by 
gain range. 
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4 



Roct -Lo c us (Int eracti ve Ex ample ) 

The following interactive example is of the open 
loop transfer function for an aircraft with an autopilot in 
the longitudinal mode described in Pole-Zero interactive 
section . 



a. Example of Graphic Response (Terminal Session) 

BEGIN EECORDIN G OF TERMINAL SESSION 
R; T = 0. 0 1/0 . 02 20:28:37 
EXECUTION BEGINS... 



OPTGRAPH 

DO YOU DESIRE TO ANALYZE: 

1. OPEN LOOP TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER OPTION NUMBER. 

1 



OPTGRAPH 

THE OPEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS 
CONTAINS : 

1. NUMBER OF STATES = 4 

2. NUMBER OF CONTROLS (INPUTS) = 1 

3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 

4. MARKOV PARAMETER = 10.0**- 6 

CLEAR SCREEN TO CONTINUE 
OPTGRAPH 

AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST , 
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AND/CE NICHOLS PLOT IS DESIRED FOR: 
INPOT # ? 

1 

OUTPUT # ? 

1 



OPTGRAPH 

OPEN LOOP TRANSFER FUNCTION 
INPUT # = 1 

CUTPUT # = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. POLE-ZERO HAP 

2. EOCT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
OPTGRAPH 

ENTER OPTION NUMBER. 

2 



ROOT-LOCUS 

OPEN LOCP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER OPTION NUMBER. 

1 
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EOOT-LOCUS 

OPEN 100F TRANSFER FUNCTION 
FIOTTIE SELECTION ANE PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTION OF TEK618) 

WHICH PLOTTER DO YOU DESIRE: 

1. TEK618 

2. VERSATEC 
CHOOSE OPTION 1 OR 2 

2 

PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 

HEIGHT = 

6 

WIDTH = 

8 

PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 

X MINIMUM = ? 

-15 

X MAXIMUM = ? 

5 

Y AXIS (IMAGINARY AXIS) 

Y MINIMUM = ? 

-5 

Y MAXIMUM = ? 

5 

PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 

X MINIMUM = -15.00 

X MAXIMUM = 5.00 

Y AXIS (IMAGINARY AXIS) 
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Y MINIMUM = -5.00 

Y MAXIMUM = 5.00 

DO YCU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES 11 OR »NC". 
n 



GAIN RANGE FOR ROOT-ICCUS PLOT 
MINIMUM GAIN = ? 



0 



500 



MAXIMUM GAIN = ? 



GAIN RANGE FOR ROOT-ICCUS PLOT 

MINIMUM GAIN = O.OOOD+OO 

MAXIMUM GAIN = 0.500D+03 

DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES” OR "NO”, 
n 



PLOT HEADING 

HOW MANY LINES OE HEADING DO YOU DESIRE ? 
{ 3 MAX ) 

ENTER NUMBER OF LINES. 

(ENTER "0" FOR NO HEADING) 

3 

A MAXIMUM OF 32 CHARACTERS PER LINE IS 

ALLOWED 

LINE 1 IS: 

aircraft with an autopilot in 
LINE 2 IS: 

the longitudinal mode 
LINE 3 IS: 
ref. 7; page 315 

PLOT HEADING 
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LINE 1 : AIRCRAFT WITH AN AUTOPILOT IN 

LINE 2 : THE LONGITUDINAL MODE 

LINE 3 : REF. 7; PAGE 315 

DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES” OR "NO", 

n 

>USIN G A PRE-ALLOCATED DATASET FOR UNIT FT17F001 
>USI NG A PRE-ALLOCATED DATASET FOR UNIT FT18F001 



A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 



ROOT-LOCUS 

OPEN LOOP TRANSFER FUNCTION 

DO YOU DESIRE TO MAKE ANY CHANGES TO: 

1. PLOTTER / PAGE SIZE 

2. GRAPH LIMITS 

3. HEADING 

4. GAIN RANGE 

5. NO ADDITIONAL CHANGES - PLOT ROOT-LOCUS 

6. NO CHANGES - EXIT ROOT-LOCUS PLOTTING 
ROUTINE 

ENTER CPTICN NUMBER. 

6 



RCOT-LOCUS 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER CPTICN NUMBER. 

3 
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OPTGRAPH 



OPEN LOOP TRANSFER FUNCTION 
INPUT # = 1 

OUTPUT # = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. POLE-ZERO MAE 

2. ROOT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPCT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
OPTGRAPH 

ENTER CPU CN NUMBER. 

7 



OPTGRAEH 

DO YOU DESIRE TO ANALYZE: 

1. OPEN LOOP TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER CPU CN NUMBER. 

4 

R; T=12. 36/15. 14 20: 22: 54 
record off 

END RECORDING OF TERMINAL SESSION 
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The preceding example follows at the end cf the 
chapter as figure 3. 2 . 

h. Example of tabular data (Terminal Session) 

The entering and exiting option menus for the 
tabular data terninal session are identical to the graphic 
response terminal session option menus, and have been elimi- 
nated from the following terminal session. 

BEGIN RECORDING OF TERMINAL SESSION 
R; T = 0. 01/0.02 20:42 :59 
EXECUTION BEGINS... 



ROOT-LOCUS 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER CPTICN NUMBER. 

2 



ROOT-LOCUS 

OPEN LOOP TRANSFER FUNCTION 

DO YOU DESIRE TABULAR OUTPUT TO GO TO: 

1. SCREEN 

2. PRINTER 

3. DISK (OPGRAPH LISTING) 

ENTER CPTICN NUMBER. 

1 



GAIN RANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN = ? 



1 

5 



MAXIMUM GAIN = ? 
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GAIN RANGE FOR ROOT-ICCOS PLOT 

MINIMUM GAIN = 0.100D+01 

MAXIMUM GAIN = 0.500D+01 

EO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES’' OR "NC". 
n 

HOW MANY POINTS DO YOU WANT TO TABULATE? 

(500 IS THE MAXIMUM) 

ENTER NUMBER OF POINTS. 

5 

ROOT-LOCUS 

OPEN LOOP TRANSFER FUNCTION 
INPUT NUMBER = 1 

CUTEUT NUMBER = 1 

TRANSFER FUNCTION (DC) GAIN = 0.1000D+01 

OPEN LOOP TF POLES AND ZEROS 



EE NOMINATOR ORDER 
POLE LOCATIONS 
REAL PART 
O.OOOOD+OO 
0. 1000D + 0 1 
-0.2000D+0 1 
-0.2000D+0 1 

NUMERATOR ORDER = 
ZERO LOCATIONS 
REAL PART 
-0.1000D+01 

GAIN = 0. 1000 

REAL PART 
0 . 7074D- 0 1 
0.8962D+00 
-0.1983D+01 



= 4 

IMAGINARY PART 
0.0000D+00 
0.0000D+00 
0.3464D+01 
-0.3464D+01 

1 

IMAGINARY PART 
0.0000D+00 

1 

IMAGINARY PART 
0.0000D+00 
0.0000D+00 
0.3441D+0 1 



E+0 



4 1 



-0.1 983D+0 1 



-0. 3441D+01 



GAIN = 0.1800 

REAL PART 
0.1 453D + 00 
0.7948D+00 
-0.1970D+01 
-0.1970D+01 

GAIN = 0.2600 

REAL PART 
0.2576D+ 00 
0 . 6552D+ 00 
-0.1 956D + 0 1 
-0.1956D+01 



1 

IMAGINARY PART 
0. 0000D+00 
0. 0000D+00 
0. 3422D+0 1 
-0. 3422D+01 

1 

IMAGINARY PART 
0.0000D+00 
0.0000D+00 
0. 3402D+0 1 
-0.3402D+01 



E+0 



E+0 



GAIN = 0.3400 

REAL PART 
0. 4426D + 00 
0 . 4426D + 00 
-0.1 943D+0 1 
-0. 1943D+01 



1 

IMAGINARY PART 
0. 1661D+00 
-0. 1661D+00 
0.3383D+01 
-0.3383D+01 



E+0 



GAIN = 



0. 4200 E+0 1 



REAL PART 
0.4285D+00 
0. 4285D+00 
-0 . 1 929D + 0 1 
-0.1 929D + 0 1 



IMAGINARY PART 
0.30 96D+00 
-0.30 96D+00 
0. 3363D+0 1 
-0. 3363D+01 



DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NO". 

Y 



ROOT-LOCUS 

OPEN LCCP TRANSFER FUNCTION 
DO YOU DESIRE TO MAKE ANY CHANGES TO: 
1. OUTPUT DEVICE 
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2. GAIN RANGE 

3. NUMBER OF POINTS CALCULATED 

4. NO ADDITIONAL CHANGES - TABULATE DATA 

5. EXIT ROOT-LOCUS TABULAR DATA ROUTINE 
ENTER CPTICN NUMBER. 

5 

ROOT-LOCUS 

OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER CPTICN NUMBER. 

3 

R; 1=1.63/3.40 20:44:54 
record off 

END RECOFDING OF TERMINAL SESSION 



G. FREQUENCY RESPONSE ANALYSIS 

The calculations for the frequency response of the Eode, 
Nyguist, and Nichols responses are identical. The only 
difference is in the graphical presentation of the calcu- 
lated data. The subroutines BODE, NYQST, and NICHOLS call 
the subroutine FREQ to calculate the transfer function 
magnitude and phase for the frequency range specified by the 
user . 

The subroutine FREQ eliminates the extraneous zeros by 
comparing the zeros to the Markov parameter as described in 
the Characteristic Equation Analysis section. The frequency 
response magnitude and phase are calculated using double 
precision complex library functions following the 
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conventions of the IBM System /360/370/ FORTRAN IV language. 
The frequency response is calculated by determining the 
product of successive numerator zeros which have been 
combined with the frequency and divided by the product of 
successive denominator poles which also have been combined 
with the frequency. The magnitude for the complex result of 
the zeros, poles, and frequency calculations is determined 
using the FORTRAN double precision library functions for 
finding the magnitude of a complex argument. The phase of 
the of the complex result is determined by separating the 
real and imaginary parts of the complex results, and finding 
the phase angle using the double precision FORTRAN library 
function arctangent. The angle is converted from radians to 
degrees and is corrected for proper quadrant by checking the 
signs of the real and imaginary parts. The normal range for 
calculating phase angles by the subroutine is from +540 
degrees to -540 degrees. This range is considered adequate 
for most application, but should these limits be exceeded 
the plot will jump with a vertical line to either +360 
degrees or -360 degrees respectively before continuing with 
the plot. 

For a graphic response the user specified frequency 
range is divided into 500 equally spaced points (on a 
logrithmic scale) for which to calculate the magnitude and 
phase. The user is given the option to select between 1 and 
500 points for tabular data. 

1 . Eode Re spon se 

The graphic response section of the subroutine BODE 
creates two plots. The first plot is the response magnitude 
in decibels versus the frequency in radians per second. The 
magnitude from the subroutine FREQ is converted to decibels 
and plotted with its corresponding frequency. The second 
plot is of the phase in degrees versus frequency in radians 
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per second. The phase from the subroutine FREQ is plotted 
with its corresponding frequency. The frequency for both 
plots is plotted on a logrithmic scale. The scale range for 
both is set automatically to include the entire range of 
data to te plotted by the subroutine MINMAX which scans the 
data for the maximum and minimum values. System information 
(tranfer function input number, output number, and (EC) 
gain) is listed in the graph legend. 

The magnitude in decibels and the phase in degrees 
is searched for "0" decibel and 180 degree crossover points 
to calculate the phase and gain margins. 

The tabular response section gives the user the 
option to select the the number of frequency points to be 
calculated between 1 and 500 for the frequency range he has 
specified. For tabular data response the user is given 
system information (numerator order, denominator order, and 
transfer function (DC) gain) for the selected transfer func- 
tion input number and output number, and the frequency 
(radians per second) , magnitude (decibels) , phase (degrees) , 
and the real and imaginary frequency response parts. 

2. Bode Respon se ( Interac tive Example ) 

The following example of compensator transfer function with 
a filter and regulator was obtained from [Ref. 8,pp. 382 - 
384 ]. 

The compensated system is represented by: 



[C ] = [ -254. 1 -19.57 ] 



CD] = co: 

[F] = | 0. 1.0 | 

1 0. -4.6 | 
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[G] =10. | 

1 0. 787 | 

[GAM ] = | 0. 1 

I 0.1 | 

[H] * [ 1. 0. ] 

[K] = 1 95.4 | 

1 4561. 1 

[Q] = £10. ] 



a. Example cf Graphic Response (Terminal Session) 

BEGIN RECORDING OF TERMINAL SESSION 

R; T=0. 0 1/0 .02 1 8:30:23 

graphics 

EXECUTION BEGINS. . . 

OPT GRAPH 

DO YOU DESIRE TO ANALYZE: 

1. OPEN LOOP TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER OPTION NUMBER. 

3 



OPTGRAPH 
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THE COMPENSATOR TRANSFER SYSTEM OBTAINED F EOM 



OPTSYS CONTAINS: 

1 . NUMBER OF STATES = 2 

2. NUMBER OF CONTROLS (INPUTS) = 1 

3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 

4. MARKOV PARAMETER = 10.0**- 6 

CLEAR SCREEN TO CONTINUE 
OPTGRAPH 

AN OPEN LCOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
AND/OR NICHOLS PLOT IS DESIRED FOR: 

INPUT # ? 

1 

OUTPUT # ? 



OPTGRAPH 

COMPENSATOR TRANSFER FUNCTION 
INPUT # = 1 

OUTPUT # = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. EOLE-ZERO MAP 

2. ROOT-LOCUS 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 

ENTER OPTION NUMBER. 
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3 



BODE PLOT 

COMPENSATOR TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

• 3. QUIT SUBROUTINE 
ENTER OPTION NUMBER. 

1 



BODE PLOT 

COMPENSATOR TRANSFER FUNCTION 
PLOTTER SELECTION AND PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTION OF TEK618) 

WHICH PLOTTER DO YOU DESIRE: 

1. TEK618 

2. VERSATEC 
CHOOSE OPTION 1 OR 2 

2 

PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 

HEIGHT = 

5.65 

WIDTH = 

8.5 

PLOTTING LIMITS FOR GRAPH 

FREQUENCY RANGE - RADIANS PER SECOND 
MINIMUM FREQUENCY = ? 

0 . 1 

MAXIMUM FREQUENCY = ? 

1000000 

PLOTTING LIMITS FOR GRAPH 

FREQUENCY RANGE - RADIANS PER SECOND 
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MINIMUM FREQUENCY = 0.100E+00 

MAXIMUM FREQUENCY = 0.100E+07 

DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES" OR "NO", 
n 

PLOT HEADING 

HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3 MAX ) 

ENTER NUMBER OF LINES. 

{ ENTER "0" FOR NO HEADING ) 

2 

A MAXIMUM OF 32 CHARACTERS PER LINE IS 
ALLOWED 
LINE 1 IS: 
filter simulation 
LINE 2 IS: 

ref. 8; pages 332 - 334 
PLOT HEADING 

LINE 1 : FILTER SIMULATION 

LINE 2 : REF. 8; PAGES 332 - 334 

DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES" OR "NO". 



A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 



GAIN MARGIN AND PHASE MARGIN 
PHASE CROSSOVER FREQUENCY IS 75.1184 RAD/SSC 
GAIN MARGIN IS -59.8906 DB 
GAIN CROSSOVER FREQUENCY IS ********** RAD/SSC 
PHASE MARGIN IS -89.9480 DEG 
CLEAR SCREEN TO CONTINUE 



BODE PLOT 

COMPENSATOR TRANSFER FUNCTION 

DO YOU DESIRE TO MAKE ANY CHANGES TO: 

1. PLOTTER / PAGE SIZE 

2. GRAPH LIMITS (FREQUENCY RANGE) 

3. HEADING 

4. NO ADDITIONAL CHANGES - PLOT BODE PLOTS 

5. NO CHANGES - EXIT BODE SUBROUTINE 
ENTER OPTION NUMBER. 

5 

BODE PLOT 

COMPENSATOR TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER OPTION NUMBER. 

3 



OPT GRAPH 

COMPENSATOR TRANSFER FUNCTION 
INPUT # = 1 

OUTPUT # = 1 

DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. POLE-ZERO MAP 

2. ROOT-LOCUS 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
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TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 

ENTER OPTION NUMBER. 

7 



OPTGRAPH 

DO YOU DESIRE TO ANALYZE: 

1. OPEN LOOP TRANSFER FUNCTION 

2. NOISE TRANSFER FUNCTION 

3. COMPENSATOR TRANSFER FUNCTION 

4. EXIT OPTGRAPH 
ENTER OPIION NUMBER. 

4 

R; T=7 . 0 1/8 .84 1 8: 39: 11 
record off 

END RECORDING OF TERMINAL SESSION 

The preceding example follows at the end of the 
chapter as figure 3.3 and figure 3.4 . 

b. Example of Tabular Data Listing File 

The terminal session for the Bode tabular data response is 
similar to the Bode graphic response terminal session. 
Appendix D contains an example of tabular output sent to a 
disk (OPGRAPH LISTING) for five points between 10 and 100 
radians/ second. 

3 . Nyg uis t Response 

The graphic response section of the subroutine NY QST 
creates a polar plot from the magnitude and phase data 
received from the subroutine FREQ. The plot is automati- 

cally scaled from a minimum radius of one by using the 
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subroutine MINMAX to determine the maximum phase and magni- 
tude values. The gain / phase margin calculations and the 
tabular data output are identical to the subroutine BODE. 

4 . N yguis t Response (Inter active E xam pi e) 

The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 
figure 3.5 . 

5 • Nic hol s Respons e 

The graphic response section of the subroutine 
NICHCI creates a single rectangular plot of magnitude in 
decibels versus phase in degrees. The magnitude from the 
subroutine FREQ is converted to decibels and plotted with 
the corresponding magnitude. The plot is automatically 
scaled using the subroutine MINMAX to determine the maximum 
phase and magnitude values. The gain / phase margin calcu- 
lations and the tabular data output are identical to the 
subroutine BODE. 

6 . N ichol s Response (I nte ractive E xa mple) 

The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 
figure 3.6 . 
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Figure 3.1 Pole-Zero Map Example 
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Figure 3,2 Root Locus Response Example. 
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Figure 3.4 Bode Response Example (Phase). 
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Figure 3.5 Nyguist Response Example. 
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Figure 3.6 Nichols Response Example 



IV. CONC LUSI ONS AND RECOMMENDA TIONS 
A. CCNC1U SION S 

An evaluation of the computational ability of OPTSYSX 
and CPTGRAPH, the program was tested using an 82 X 82 
matrix, provided by NASA-Edwar ds, of longitudinal equations 
for the X-29A experimental forward-swept wing fighter 
aircraft prototype. The OPTSYSX program array had to be 
redimensioned and a 2-megabyte virtual machine size was 
required for a system of this size. 

The open loop eigenvalues calculated by the OPTSYSX 
program compared favorably with the Eigensystem data 
supplied by NASA-Edwards, but the OPTSYSX program failed in 
calculating the open loop transfer function zeros. 
Extraneous zeros (calculated by OPTSYSX) were of the same 
approximate magnitude as the open loop system zeros. This 
prevented the identification and elimination of the extra- 
neous zeros by comparison with the Markov parameter. The 
open loop zero identification was further complicated by not 
having available open loop transfer function zero informa- 
tion from NASA-Edwards to compare the OPTSYSX results to. 

When requested for additional information NASA-Edwards 
provided a computer tape with a revised 98 X 98 matrix of 
longitudinal motion equations for the X-29A fighter 
aircraft. For this revised system, the OPTSYSX program 
requires a virtual machine size in excess of 2-megabyte. 
Time constraints have prevented running the OPTSYSX program 
with the revised system. 

The OPTGRAPH program was tested using the OPTSYSX 
program’s open loop Eigenvalue and transfer function output 
data for the X-29A 82 X 82 system of longitudinal motion 
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equations. OPTGF.APH provided excellent quality Bode plots 
for all of the open loop transfer function inputs and 
outputs and did not exhibit any computational difficulties 
as indicated by the lack of error messages at the end of the 
test. While the test results are encouraging as to the 
ability of the OPTGRAPH program to assist in the analysis of 
large order systems, the actual computational accuracy of 
the OPTGRAPH program can not be fully verified until accu- 
rate transfer function information can be obtained from the 
OPTSYSX program. 

B. RECOMMENDATIONS 

Eased on the results of this Thesis the following areas 
emerged as possible areas for further research and study. 

1 . Root Fi n ding Routin es 

The OPTSYSX program has two similar root finding 
subroutines (HQR and HQR2) , to find the transfer function 
Eigensystem (poles) and zeros. The order of the transfer 
function numerator is calculated separately in the subrou- 
tine ZEROS. 

The inclusion of the open loop transfer function 
eigensystem (pole locations) , zero locations, transfer func- 
tion numerator orders, and transfer function gains, for the 
revised (98 X 98) longitudinal mode system by NASA-Edwards 
should simplify the locating the problem in calculating the 
transfer function zeros. 

The open loop system data from NASA-Edwards also 
included Bode magnitude and phase plots, and gain and phase 
margin information, which can be used to evaluate more 
conclusively the OPTGRAPH program’s ability to help evaluate 
a large order system. 
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2 • Fr og ra m Mem ory Requi rements 

The OPTSYSX program requires in excess of 2-megabyte 
of virtual machine memory when configured for large matrix 
operations (98 X 98) . Virtual machines with this memory 
capacity are not normally available to a user. The memory 
usage for matrix storage and manipulation is a possible area 
for the reduction in program memory requirement size. The 
OPTSYSX program retains most of its computational arrays in 
memory. This method simplifies programming but is extremely 
inefficient in its use of virtual machine memory. 

3 ♦ F urth er Modi f ica tion s 

The area of modern digital controls should be further inves- 
tigated. The computational abilities of the OPTGRAPH 
program make it readily adaptable, with minor modifications, 
for assisting in the analysis of transfer functions in the 
digital domains. 
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APPENDIX A 
OPTSYS EXEC LISTING 



5TRACE OFF 

************ ****************** ***************** ******** **** 

* THE OPTSYS EXEC 

* CONTROLS THE OPTSYSX, OPTCALC AND OPTPLOT 

* TO DETERMINE THE TIME RESPONSE OF A SYSTEM. 

* BY H. A. DIEL 

* VERSION 1.0 16 JULY 1984 

* MODIFIED BY M. H. LAPTAS 

CLRSCEN 

SBEGTYPE -ENDFOUR 



**** * ****** *** * *** * ****************************** ******** ** 

OUTPUT FROM OPTSYSX HILL COME TO THE SCREEN IF YOU 
WISH DISK FILE 

ENTER D ANY OTHER INPUT YIELD SCREEN 
****** ********* ******************************************** 

-ENDFOUR 

BREAD VARS BANS 

SIF .SANS NE .D SGOTC -ONE 

FILEDEF 06 DISK OUTPUTX LISTING A1 

CLRSCRN 

BBEGTYPE -ENDSIX 
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*^********* ************************************ ************ 



OUTPOT WILL GO TO DISK FILE •OUTPUTX LISTING A1 f 
*********************************************************** 

-ENDSIX 

CP SLEEP 3 SEC 

-ONE 

CLRSCRN 

CP LINK 0039P 191 251 HR 
ACC 251 E 

*********************************************************** 
* * 

* ALLOW THE USE CE AN OLD "OPTMAT DATA A1» * 

* * 
****************************************** ***************** 
RENAME OPTMAT DATA A1 OPTSYS DATA A1 

SIF SRC NE 0 SGOTO -EIRST 

RENAME OPTSYS DATA A1 OPTMAT DATA A1 

FINDSTAK OPTMAT DATA A1 C01 0 LIM002 ALL G ROUP 1 

SREAD VARS SA 1 SA2 

SIF . S A2 NE .1 SGOTO -START 

CLRSCEN 

SBEGTYPE -ENDTWO 

*********************************************************** 

YOU HAVE A DATA FILE NAMED ’OPTMAT DATA' ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 
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IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 
THIS RUN , THE CPTSYS PROGRAM WILL READ IN THE 



DESIRED DATA AT THE APPROPRIATE TIME/ 

IF YOU TYPE (Y) ES. 

ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 

£******$*$$ ** ****$***$**:!=**:!«:!«:<«*********:$::!«:!=$ * # * * * * * * * ** 

-ENDTWO 

SREAD VARS SANS 

SIF .SANS EQ .Y SGOTO -START 

-FIRST 



* * 

* ERASE THE OLD "OPTMAT DATA A1 " DATA FILE * 

* PLACE ”000 0" IN THE NEW "OPTMAT DATA FILE" * 

* TO ACT AS A FLAG FOR OPTSYSX AND OPTCALC * 

* * 



ERASE OPTMAT DATA A1 
SSTACK 000 0 

FILESTCK OPTMAT DATA A1 F 80 1 

-START 

CLRSCFN 

SBEGTYPE -ENDONE 

THE OPTSYS EXEC CONTROLS FOUR PROGRAMS: 

1 OPTSYSX FORTRAN (SYSTEM ANALYSIS) 

GENERATES OPTMAT AND OPTGRAPH DATA SETS 
- RERUN USES OPTMAT 

2 OPTCALC FORTRAN (CALCULATE TIME RESPONSE) 

GENERATES OPTPLOT DATA SET USES OPTMAT FROM 
OPTSYSX 

3 OPTGRAPH FORTRAN (POLE-ZERO/ ROOT-LOCUS , BODE, 
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NYQUIST, ETC) 

USES 0PTGR01 -OPTGRNO -OPTGRCM FROM OPTSYSX 
4 EXIT 

OPTPLOT FORTRAN IS THE FOURTH PROGRAM 
USES OPTPLOT DATA SET FROM OPTCALC 
TEE SIZE OF THE DATA SETS VARY WITH THE SYSTEM 
ORDER, AND CAN USE ABOUT 20% OF THE USERS DISK SPACE. 
THEREFORE ENSURE THAT SUFFICIENT DISK SPACE' IS 
AVAILABLE. 

IF DATA SET IS ALREADY AVAILABLE YOU MAY RUN ANY 
PROGRAM 



- ENTER 1 



ANY OTHER INPUT 
RETURNS TO MENU 



-ENDONE 

SREAD VARS SANS 

SIF .SANS EQ .1 SGOTO -OPTSYS 

SIF .SANS EQ .2 SGOTO -OPTCALC 1 

SIF .SANS EQ .3 SGOTO -OPGRAPH 1 

SIF .SANS EQ .4 SEXIT SRC 

SGOTO -START 

-OPTSYS 

FILEDEF 8 DISK OPTPLC1 DATA Al (PERM 
FILEDEF 9 DISK OPTMAT DATA Al (PERM 
FILEDEF 10 DISK OPTGECL DATA Al (PERM 

FILEDEF 01 DISK OPTGRNO DATA Al (PERM 

FILEDEF 04 DISK OPTGRCM DATA Al (PERM 

GLOBAL TXTLIB VFORTLIE IMSLDP NONIMSL 

STYPE LOADING OPTSYS GENERAL STATE VARIABLE ANALYSIS 

PROGRAM 

STYPE OPTMAT DATA Al INPUT DATA SET ON RERUNS 

OPTSYSX 

CLRSCEN 

STYPE... DATA SETS .. OPTMAT AND .. OPGRAPH CREATED 
CP SLEEP 5 SEC 
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6G0T0 -START 
-0PTCALC1 



* * 

* CHECK FOE DATA IN THE FILE "OPTMAT DATA " * 

* EEFORE LOADING OPTCALC * 

* * 



sit####****:*!**###*##**#*###*#*;}:*####**######:***#*#*#*### jfc## 

FINDSTAK OPTMAT DATA Al C01 0 LIM002 ALL GROUPl 

DREAD VARS SA 1 SA2 

SIF . £A2 EQ .1 5 GOTO -OPTCALC 

5TYPE PROPER DATA FILE IS NOT AVAILABLE FOR OPTCALC 
CP SLEEP 5 SEC 
5GOTO -START 
-OPTCALC 

5TYPE OPTCALC TIME RESPONSE PROGRAM 

FILEDEF 8 DISK OPTPLCT DATA Al (PERM 
FILEDEF 9 DISK OPTMAT DATA Al (PERM 
GLOBAL TXTLIB VFORTLIE IMSLDP NONIMSL 
OPTCALC 

STYPE .. OP T PLOT EATA Al CREATED 

C? SLEEP 5 SEC 

CLRSCEN 

SBEGTYPE -ENDNINE 






IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE * Y* TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 
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*********************************************************** 



-ENDNINE 

SREAD VARS SANS 

SIF .SANS EQ . Y SGOTC -START 

STYPE PLOTS NEXT 

CP SLEEP 4 SEC 

GLOBAL TXTLIB DI SLIBVS 92DISLIB 92INTLIB VFORTLIB GRFLIE 
NONIMSL CMSLIB 

FILEDEF 8 DISK OPTPLOT DATA A1 (PERM 

FILEDEF 11 DISK LBLAKK TMP A3 (RECFM F LRECL 2400 BLKSIZE 
2400 ELKSIZE 2400 XTENT 600 
FILEDEF 12 DISK ISSCCMAP MAPDTA * (RECFM F LRECL 400 
BLKSIZE 400 

FILEDEF 13 DISK TABLET TMP A3 (RECFM VS LRECL 208 BLKSIZE 
208 XTENT 1000 

FILEDEF 14 DISK LBLANK MAPDTA SLBMODE (RECFM VS LRECL 608 
BLKSIZE 608 

FILEDEF 17 DISK DISSPIA SYSOT1 A3 (RECFM FB LRECL 2000 
BLCCK 2000 XTENT 500 

FILEDEF 18 DISK DISSPIA METAFILE A4 (RECFM VBS LRECL 19065 
BLCCK 19069 
OPTPLOT 
SGOTO -META 
-OPGRAPH1 

*** *************************************************** *** ** 



* * 

* CHECK FOE DATA IN THE OPTGRAPH DATA SETS " * 

* EEFORE LOADING OPGRAPH * 

* * 



*** * ************************************** ***** ** * **** * ** ** 
*********************************************************** 
RENAME OPTGROL DATA A1 OPTSYS DATA A1 
SIF SEC NE 0 SGOTO -CFTNOISE 
RENAME OPTSYS DATA A1 OPTGROL DATA A1 
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FINDSTAK OPTGROL DATA Al C01 0 LIM002 ALL GROUP 1 
GSEAD VARS GAl &A2 S A3 GA4 &A5 
GIF . SA2 NE .0 &GOTO -SIXTEEN 
-OPTNCISE 

RENAME OPTGRNO DATA Al OPTSYS DATA Al 

GIF &RC NE 0 &GOTO -CPTCOMP 

RENAME OPTSYS DATA Al OPTGRNO DATA Al 

EINDSTAK OPTGRNO DATA Al C01 0 LIM002 ALL GROUP1 

GREAD VARS GA1 &A2 &A3 GA4 &A5 

GIF . &A2 NE .0 &GOTO -SIXTEEN 

-OPTCCMP 

RENAME OPTGRCM DATA Al OPTSYS DATA Al 

GIF &RC NE 0 &GOTO - NCDATA 

RENAME OFTSYS DATA Al OPTGRCM DATA Al 

EINDSTAK OPTGRCM DATA Al C0 1 0 LIM002 ALL GROUP 1 

SREAD VARS SA 1 SA2 SA3 SA4 &A5 

GIF . &A2 NE .0 SGOTO -SIXTEEN 

-NODATA 

SPRINT PROPER DATA FILE IS NOT AVAILABLE FOR OPGRAPH 
C? SLEEP 5 SEC 
SGOTO -START 
-SIXTEEN 

STYPE ...OPTGRAPH CLASSICAL ANALYSIS OF OPTSYS OUTPUT 

SET LDRTELS 10 

GLOBAL TXTLIB DISLIBVS 92DISLIB 92INTLIB V FORTLIB GRFLIB 
NONIMSL CMSLIB 
SET LDRTBLS 10 

FILEDEF 01 DISK OPTGRNO DATA Al (PERM 
FILEDEF 02 DISK OPGRAPH LISTING Al 
FILEDEF 03 PRINTER (RECFM FA BLKSIZE 133 PERM 
FILEDEF 04 DISK OPTGRCM DATA Al (PERM 
FILEDEF 06 TERM (RECFM FA BLKSIZE 133 
FILEDEF 05 TERM (RECFM FA BLKSIZE 80 
FILEDEF 10 DISK OPTGECL DATA Al (PERM 

FILEDEF 11 DISK LBLANK TMP A3 (RECFM F LRECL 2400 BLKSIZE 



68 



2400 XTENT 600 



FILEDEF 12 DISK ISSCCMAP MAPDTA * (RECFM F LRECL 400 
BLKSIZE 400 

FI1EDEF 13 DISK TABLEI TMP A3 (RECFM VS LRECL 208 BLKSIZE 
208 XTENT 1000 

FILEDEF 14 DISK LBLANK MAPDTA SLBMODE (RECFM VS LRECL 608 
BLKSIZE 608 

FILEDEF 17 DISK DISSPLA SYSUT1 A3 (RECFM FB LRECL 2000 
BLOCK 2000 XTENT 500 

FILEDEF 18 DISK DISSPLA METAFILE A4 (RECFM VBS LRECL 19065 
BLOCK 19069 
CPTGRAPH 
-META 



* * 

* CHECK FOR FILE "DISSPLA METAFILE A4" ON * 

* THE USER’S DISK BEFORE GOING TO DISSPOP * 

* * 



RENAME DISSPLA METAFILE A4 OPTSYS METAFILE A4 
SIF SRETCODE NE 0 SGCTO -START 

RENAME OPTSYS METAFILE A4 DISSPLA METAFILE A4 

-EIGHT 

CLRSCRN 

SBEGTYPE -ENDTEN 



DC YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
THE DISSPLA METAFILE THAT YOU JUST CREATED? 



(Y OR N) 



****************************************** ****************** Z* 

-ENDTEN 

SREAD VARS SANS 

SIF .SANS EQ .Y SGOTC -NINTH 

SIF .SANS EQ . N SGOTC -START 

SGOTO -EIGHT 

-NINTH 

EXEC DISSPOP VRSTEC 
CLRSCRN 

SBSGTYPE -ENDTWELVE 



YOUR GRAPH (S) CAN EE PICKED UP AT THE COMPUTER CENTER. 

THE GRAPH (S) Will BE ADDRESSED TO "POP (USER ID)' 1 . 

****************************************************** ## 

-ENDTWELVE 
CP SIEEP 5 SEC 
SGOTO -START 
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APPEN DIX B 

OPTSYSX PROGRAM LISTING 

This portion of the thesis contains the OPTSYSX FORTRAN 
program (88 pages) . 
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OPTGRAPH PROGRAM LISTING 

This portion of the thesis contains the OPTGRAPH FORTRAN 
program (54 pages) . 
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